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By cooling of equilibrium lattice fields at finite temperature in SU(2) gauge theory it has been 
shown that topological objects (calorons) observed on the lattice in the confined phase possess a 
dyonic substructure which becomes visible under certain circumstances. Here we show that with 
the increasing temporal lattice extent the distribution in the caloron parameter space is modified 
such that the calorons appear non-dissociated into constituent dyons. Still the calorons have non- 
trivial holonomy which is demonstrated by the Polyakov line behaviour for these configurations. At 
vanishing temperature (on a symmetric lattice) topological lumps obtained by cooling show rota- 
tional symmetry in AD for the action density, but a characteristic internal double peak structure of 
Polyakov lines with respect to all (temporal and spatial) directions. 

PACS numbers: ll.10.Wx, 12.38.Lg, 14.80.Hv 



I. INTRODUCTION 

At high temperatures, near but below the deconfmc- 
ment temperature, classical solutions of Yang-Mills equa- 
tions with nontrivial holonomy (KvB calorons 0,0) are 
seen on the lattice for SU{2) gauge theory to be fre- 
quently dissociated into dyons 0, 0, IS 01 This means 
that the distance between the dyons forming a caloron, 
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is larger than the size of dyons (which is b/n for the 
case of maximally nontrivial holonomy). Here p is the 
instanton size parameter and b is the temporal periodicity 
interval. 

These observations have been made by the use of cool- 
ing. Therefore, the question arises why such a substruc- 
ture has not been observed by previous authors who have 
used the cooling method. In this paper we will argue that 
there is only a certain window of temperature or space- 
time asymmetry where it can be revealed by this method. 

The interest in the existence of caloron constituents has 
increased since it has been demonstrated in Ref. 7] that a 
constituent substructure very reminiscent of the caloron 
solutions can also be identified without cooling, above 
and below the phase transition. This can be achieved 
by using the localization properties of the fermionic zero 
modes of a suitably chirally-improved Dirac operator. 
The similarity with the properties of a caloron solution 
is strikingly realized for a certain fraction of configura- 
tions with topological charge Q = ±1, where the sin- 
gle zero mode is seen to change its localization when 
the periodicity of fermionic boundary conditions becomes 
modified. A systematic study |jj of the typical pattern 
of localization and derealization followed by jumps of 
the zero mode has revealed that this pattern depends 



on the timelike holonomy exactly in a caloron-like way. 
Whereas the topological density has a much more com- 
plicated structure, the positions where the zero mode is 
pinned-down actually show the signatures expected for 
caloron constituents 8J: they are local maxima of the 
topological density q{x) with a sign as required by the 
chirality of the mode. This suggests that (dissociated 
or non-dissociated) calorons might really form the semi- 
classical background of the gauge field near the phase 
transition. 

Coming back to SU(2) calorons with their two con- 
stituents, it seems that the quantities d and p appearing 
in eq. Q are impossible to be measured simultaneously: 
when d is seen by observation of separate dyon positions 
no instanton-like profile (of topological density) is ob- 
served which could be used to define p. When d goes to 
zero (d << b) an instanton size parameter p can be mea- 
sured by comparison with the instanton's action density 
profile, but then no dyons are seen as separate objects. 
More precisely: the parameter d cannot be measured for 
all caloron configurations as the distance between con- 
stituents as long as only the action or topological charge 
densities are available as local observables to describe 
them. 

The time periodicity parameter b defines the temper- 
ature T: b = 1/T. In order to demonstrate how the re- 
combination of constituents depends on the temperature 
we can change the temporal extent b of the lattice. 

We employ the standard relaxation cooling technique 
using the Wilson lattice action and concentrate on the 
investigation of lowest-action field configurations with 
Q = ±1. Besides of the fact that solutions with \Q\ = 1 
are not existing on a torus in a mathematically strict 
sense [i| , also the cooling method has a limited relevance. 
It cannot be used for revealing the full topological struc- 
ture of the QCD vacuum even if the latter is semi-classical 
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to a certain extent. In particular, it considerably weak- 
ens the chiral condensate compared with its value for 
equilibrium fields ^j- Moreover, it is well-known that 
the Wilson action depends on the instanton (or caloron) 
scale-size. Therefore, its minimization shrinks the local- 
ized solutions until they disappear 'through the lattice 
meshs'. Thus, small excitations will be lost first. In the 
literature - besides fermionic methods - there are bet- 
ter techniques like improved cooling, smoothing, smear- 
ing etc. allowing to decipher the topological long-range 
structure of the gauge fields. But this task is not our 
concern in this paper. Here we ask the more modest 
question, what kinds of simple classical solutions can be 
found from equilibrium gauge fields at different temper- 
atures by successively minimizing the action. 

This question remains interesting for those who want 
to build or apply semi-classical-like approximations, i.e 
models like the instanton gas or liquid model 0] be- 
ing successful in many phcnomcnological applications for 
which chiral symmetry breaking (and not quark confine- 
ment) is playing the major role |l2(. Until now most of 
the cooling (or smoothing) results obtained at T < Tc 
have been interpreted in terms of BPST instanton |l3| 
and Harrington-Shepard (HS) caloron 0] solutions, re- 
spectively, all exhibiting trivial asymptotic holonomy. 
Here we would like to convince the reader that the 
simplest caloron or instanton solutions seen after cool- 
ing have typically non-trivial holonomy, irrespectively of 
their possible dissociation into dyons pairs. Therefore, 
they cannot be correctly interpreted as BPST or HS so- 
lutions. 

The paper is organized as follows. In section ^ we 
draw the attention to the static nature of configurations 
near the deconfining transition. Section lllll presents our 
results on semiclassical configurations at finite tempera- 
tures, pointing out the loss of "staticity" and the increas- 
ing importance of the Polyakov loop for detecting the 
non-trivial substructure at lower temperature. In section 
II VI we extend the cooling studies to the symmetric torus. 
We emphasize that the KvB caloron soultions were con- 
structed in an infinite spatial volume. When used at 
finite volume considerable deviations are to be expected 
when the spatial box size is no longer large compared to 
b, since the typical size of the constituents is b/n. Fi- 
nally, in section we discuss the consequences, also in 
the perspective of a twin paper by Gattringer et ai. 



II. NON-STATICITY AND SEPARATION INTO 
CONSTITUENTS 

It turns out that the possibility to observe the dyonic 
constituents of a KvB caloron as lumps of action depends 
on (Y) 2 - In SU(2) LGT at = 4,/gl = 2.2 on a lattice 
16 3 x 4 the parameter p is concentrated near the value 
p « 2.5a (a is the lattice spacing) Q. With b — 4a, 



This means that dyons are well separated. The value of 
p in [J] was determined by fitting the lattice caloron with 
the analytic KvB caloron, and formula JQ) was used. On 
the lattice 16 3 x 6 (with b = 6a) and at the same j3 = 2.2 
(i.e. at a temperature 1.5 times lower) the parameter 
(^) 2 would be of the order 0(1). Then, from this sim- 
ple arithmetics, one would expect that calorons are not 
dissociated into dyonic lumps anymore. 

The possibility to measure the distance between dyons 
inside a caloron just by detecting the peaks of the action 
density on the lattice is given only in the case of well- 
separated objects. Indirectly this distance can be mea- 
sured by measuring a quantity that can be called non- 
staticity [2(j. Unfortunately, cooling yields metastable 
plateaus only for temperatures below the deconfining 
temperature. On the other hand, this allows us to re- 
strict ourselves in the following to maximally non-trivial 
holonomy because the average Polyakov line vanishes. To 
describe for this special case the relation between dis- 
tance d and non-staticity, we have considered analytical 
caloron solutions in continuous space-time. We divided 
the time interval b into N t time slices and expressed the 
action in the i-th time slice, Si = Y]^ ss,i, in terms of the 
local action density sj^. The non-staticity St is defined 
as 
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Obviously, this definition depends on the number of time 
slices Nt such that one would have 5t — > for N t — > oo. In 
order to get an asymptotically Nt independent quantity 
we will modify the definition of 5t as follows 



Ei='i T,s\ s s,i+i ~ s x,i\ N t 
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The factor 1/4 has been chosen such that for a lattice 
16 3 x4, where all the simulations were initially done using 
(J2J), the two definitions © and agree. 

With the definition © at hands we can calculate the 
non-staticity of an analytic KvB caloron. The non- 
staticity depends on the holonomy and on the distance 
between the constituents inside the caloron. For max- 
imally nontrivial holonomy (which coincides with the 
average holonomy in the confinement phase) the con- 
stituent dyons have equal mass. For this simplified 
case we have determined a bifurcation value of the non- 
staticity, b~l — 0.27, choosing the distance between the 
constituents such that the two lumps of action density 
(dyons) merge into one lump of action density. This sin- 
gle lump is what we call a "recombined caloron". As 
already mentioned, the non-staticity uses the discrete- 
ness of a lattice configuration. The value S^ is obtained 
by inserting the analytic form of the action density of a 
continuum KvB caloron solution [lj, calculated exactly 
at the point of bifurcation, into the evaluation of © us- 
ing a grid with lattice spacing a = b/N t . In order to see 
how well-defined at finite Nt this "bifurcation value" 5 t * 
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FIG. 1: The dependence of non-staticity , denned at the 
bifurcation point for analytic KvB caloron solutions, on the 
number of time slices N t ■ The zigzag form of the curve can be 
explained by the qualitatively different arrangement of even 
and odd time points where the continuum KvB caloron solu- 
tion has to be calculated in order to evaluate iffi. 



can be, we have evaluated it for various (even and odd) 
N t > 4. The function 8;(N t ) that goes to <5 t * = 0.27 for 
Nt — > oo is presented in Fig. 2] 

If the non-staticity is lower than the bifurcation value 
8t two dyons can be distinguished by the two maxima of 
the action density. If the non-staticity is bigger than 5^ 
two dyons appear recombined into a caloron with only 
one action density maximum. 

III. RECOMBINATION OF DYONS INTO 
CALORONS WITH LOWERING TEMPERATURE 

The main results of this paper are obtained from en- 
sembles of SU(2) gauge field configurations created by 
heat bath Monte Carlo at j3 = 2.2 with respect to the 
Wilson action Sw on lattices 16 3 x Nt with N t = 4, 
N t — 5 and N t = 6. Each ensemble consisted of about 
8000 independent configurations . Cooling of these con- 
figurations was performed using the fastest possible re- 
laxation with respect to the Wilson action. For this 
method each link U Xt(l is immediately replaced by the 
projection to 5(7(2) of the staples around it, U xfl . 

The cooled configurations studied in this paper have 
been identified when the cooling history finally has ar- 
rived at a quasi-stable plateau on the level of a single 
instanton action Sw ~ Sjnst = 2tt 2 /3. Lateron, we de- 
scribe also some results concerning higher plateaux of 
action, e.g. Sw ~ 3Si ns t which are passed earlier in 
the cooling process. More precisely, on such plateaux we 
stopped cooling always at local minima of the violation 
of the lattice equations of motion Here, the violation 
is defined as 

/ 1 \l/2 

v = ( 2 tr [Ux » ~ [Ux » ~ ^ ) ■ ( 4 ) 



In addition to that, the following conditions have been 
imposed for the automatized selection of the classical so- 
lutions: 

• the action fits into the window 0.5 < Sw / 'Si ns t < 
1.25, 

• the decrease of action has slowed down to 
\AS W \/S inst < 0.05 and 

• the violation of the equations of motion should be 
sufficiently weak, V < 25. 

In most of the cases for which the first and second con- 
ditions hold we find a small minimum violation V < 20. 
In a few cases there were also local minima with much 
larger ^-values, which hardly could be interpreted as 
extended solutions of the lattice equations of motion. 
Thus, the third condition was implied in order to select 
"good" solutions. We have varied our stopping crite- 
ria |2H an( A seen that the details of the ensemble of the 
'frozen' objects slightly change. Our main observations 
reported below do not depend on the details of the pro- 
cedure. The efficiency of the conditions was such that 
80 % (in the case N t = 4), 60 % (N t = 5) and 55 % 
{Nt — 6) of the equilibrium configurations ended up in a 
cooled configuration at the one-instanton action plateau. 
These cooled configurations include events with topolog- 
ical charge Q = ±1, i.e. real calorons, as well as con- 
figurations with topological charge Q — which in our 
previous paper had been identified as static dyon-anti- 
dyon pairs (DD). The latter constitute 18 % (N t = 4), 
14 % (N t = 5) and 8 % (N t = 6), respectively, of all 
cooled configurations at the one-instanton level. They 
have been discarded from the considerations in this pa- 
per. 

We remind the reader that the recombination thresh- 
old 8t — 0.27, strictly speaking, reflects the recombina- 
tion for maximally non-trivial holonomy only, i.e. with 
an asymptotic value of the Polyakov line L as = 0. If one 
performs cooling without special restrictions concerning 
the holonomy, there is no guarantee that the asymp- 
totic holonomy of the caloron configurations still coin- 
cides with the average Polyakov line of equilibrium con- 
figurations in confinement. For the purpose of defining 
an asymptotic holonomy L as for each cooled configura- 
tion, we have determined the average of over a 3D 
subvolume where the local 3-D action density sg is low, 
for definiteness sg < 0.0001. In Fig. |3we present the dis- 
tribution of cooled configurations over L as as a histogram 
(with bin size 0.1) for the three cases Nt — 4, Nt = 5 and 
Nt = 6. In the legend we show the respective volume 
fraction (F « 0.15) of the three- volume over which the 
"asymptotic" value L as is defined as an average, i.e. far 
from the lumps of action and topological charge. 

As explained above, the non-staticity St is a measure 
which describes the distance from a perfectly (Euclidean) 
time independent configuration. In other words, the dis- 
tributions of non-staticity of caloron events obtained by 
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FIG. 2: The distribution of holonomy L as for the three sam- 
ples of cooled configurations corresponding to three temper- 
atures below deconfinement at Nt = 4, Nt = 5 and Nt = 6. 



cooling of equilibrium lattice configurations can be con- 
sidered as a substitute for the distribution in dyon dis- 
tances d. This quantity can be directly measured for 
cooled lattice gauge field configurations. We show in 
Fig. 01 the St distributions for all our cooling products 
obtained at /3 = 2.2 on 16 3 x Nt lattices. 

In an attempt to make a fair comparison with calorons 
with non-trivial holonomy and to correct for the possible 
evolution of the asymptotic holonomy away from L as = 
during the cooling process, we defined a subsample by the 
requirement |L as | < 1/6. One can see that the cut with 
respect to the asymptotic holonomy selects cooled config- 
urations from the flat central part of the histogram shown 
in Fig. |2 On the other hand, we notice that a consid- 
erable fraction of cooled configurations has developed an 
asymptotic holonomy \L as \ > 1/6. 

In Fig. |2| (a) we show the probability distribution 
over 6t for cooled configurations with an action at the 
one-instanton plateau without the cut according to the 
asymptotic holonomy |£ os |. One can see that a relatively 
high fraction of configurations, obtained from the Monte 
Carlo equilibrium with N t = 4, has S t < S% = 0.27. This 
means that they would be identifiable as consisting of 
two constituents by looking for the 31? action density on 
the lattice. For Nt = 5 it is only a minority of cooled 
configurations which falls below the threshold St = 0.27. 
No static (according to the non-staticity criterum) con- 
figurations have been found among cooling products at 
N t = 6. 

We have repeated the same analysis after applying the 
cut with respect to the asymptotic holonomy, |L as | < 
1/6. Then we get modified histograms in St for the three 
temperatures. This is shown in Fig. [3] (b). The his- 
tograms got more pronounced peaks in St which are po- 
sitioned around 0.125, exactly around S% = 0.27, and 



around 0.5 for Nt — 4, Nt — 5 and N t = 6, respectively. 

There are other criteria which could be used to char- 
acterize a more or less static configuration, for exam- 
ple, the presence of static Abelian monopoles emerging 
in the maximal Abelian projection. This criterium is not 
identical with the separation set by 5%. Thus, another 
subsample can be defined by the property that a pair of 
static Abelian monopoles has been found after fixing the 
cooled configuration to the maximal Abelian gauge and 
doing the Abelian projection. This subsample can be an- 
alyzed with respect to the three-dimensional distance R 
between the Abelian monopoles. 

In Fig. 01(a) we show the histogram over St of all cooled 
configurations obtained from the Monte Carlo ensemble 
at Nt = 4 together with the histogram of those which 
cxplicitely exhibit the dyon-dyon structure in terms of 
Abelian monopoles. One can see that practically all 
cooled configurations below <5jf possess this structure, but 
above S% the fraction rapidly goes to zero. We show the 
same for Nt = 6, i.e. at lower temperature, in Fig. E](b). 
In this case, at the peak value around St = 0.5 only 20 % 
of the solutions still are characterized by a static dyon- 
dyon pair, whereas at higher non-staticity this is never 
the case. In these two distributions no cut with respect 
to the asymptotic holonomy has been applied. 

Up to now we have two tentative definitions of the po- 
sition of the constituents: one is the position of the two 
static Abelian monopoles (in MAG, as long as they they 
are static) and the other is defined by the maxima of 
the 3D action density. [2^ In the case of well-separable 
maxima of the latter these maxima fall close to the posi- 
tions of the MAG monopoles. In the other limiting cases 
of recombined maxima or imbalanced maxima (this cor- 
responds to an asymptotic holonomy far from zero) of 
the 3D action density at least the absolute maximum of 
action density is still easy to find. It is either the sin- 
gle maximum or it takes the role of one of the maxima. 
For the case of Nt = 4, i.e. the temperature just be- 
low deconfinement, we show in Fig. [3] (a) the histogram 
with respect to the local Polyakov loop at the two sorts 
of 3D constituent points, the loci of static monopoles or 
the maximum of action density. In both definitions the 
histogram peaks near to Lg = ±1. The peak is how- 
ever more pronounced for the monopole locations, less 
pronounced for the maxima of the 3D action density. 

The relation between the non-staticity St and the dis- 
tance R between the static dyonic constituents emerging 
(in MAG as Abelian monopole) given in lattice units, is 
presented in Fig. [3] (b) for the higher temperature, near 
the deconfinement temperature (N t = 4). For this tem- 
perature such dyon-dyon events are clearly distinguish- 
able among the cooled configurations as long as the non- 
staticity St < 0.6. For extremely low non-staticity St (left 
from the peaks in Fig. 01(a) and (b)) we find an average 
distance R w 7, whereas near 8% the average distance 
is R m 4. For higher St, the part of solutions which still 
possesses a clear dyon-dyon structure in terms of Abelian 
monopoles, has them localized at distances R between 
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FIG. 3: The distribution of non-staticity St after cooling as histograms (with bin width 0.05), for three temperatures below 
deconfinement at N t = 4, N t = 5 and N t = 6 : (a) without any cut, (b) when a cut \L as \ < 1/6 is applied. The thick vertical 
line marks the non-staticity 51 where the caloron recombines. 
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FIG. 4: The distribution of non-staticity St after cooling compared with the subsample of configurations which have a static 
dyon-dyon (DD) structure exhibited by monopoles after Abelian projection: (a) for the higher temperature with N t = 4, (b) 
for the lower temperature with Nt = 6. The thick vertical line marks the non-staticity 51 where the caloron recombines. 



one and two lattice spacings. 

In order to represent how the Polyakov line behaves in- 
side a lump of action, we have chosen the absolute max- 
imum of the 3-D action density (denoted as central point 
x ) and have explored the Polyakov loop in its neigh- 
borhood. For this purpose, we have denned a locally 
summcd-up Polyakov line L to t (summed over the central 
point xq and its six nearest neighbors a?j (i = 1, .., 6): 

6 

Ltot = ^2 L s . (5) 

i=0 



and a kind of Polyakov-line dipole moment over the same 
set of 3D lattice points with respect to the central point: 

6 
i=l 

The absolute value \L tot \ of the first quantity tests the 
amount of local coherence of the Polyakov line. The 
absolute value |Mtot| of the second quantity tests the 
amount of presence of opposite-sign Polyakov lines rep- 
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FIG. 5: The dyon-dyon structure after Abelian projection : (a) distribution of the local Polyakov loop at the position of the 
static Abelian monopoles compared to the distribution at the maxima of action density, (b) average distance between the static 
Abelian monopoles vs. non-staticity S t assigned to the cooled configuration. The thick vertical line marks the non-staticity S* 
where the caloron recombines. 



resenting eventually two different constituents inside the 
same lump of action. Fig. 0(a) shows how the \L tot \ (i.e. 
the locally summed-up Polyakov line) changes with St in 
different bins of width 0.1. For the temperature nearest 
to the transition, at N t — 4, we see that \L to t \ falls from 
4.0 to 1.0 at St > 0.5. We interprete this such that in 
the region, where constituents can be well separated ac- 
cording to the action density (at small S t ), they are char- 
acterized by a relatively smooth change of the Polyakov 
line inside. In the region of large S t where they are not 
separable according to the action density, the Polyakov 
line changes rapidly in the neighbourhood of the absolute 
maximum of action density. For the lower temperatures, 
N t — 5 and Nt = 6, the relationship between these prop- 
erties of an action cluster and St is the same. The dif- 
ference is that separable lumps of action (those with low 
St) become very rare. Fig. El (b) shows how the "dipole 
moment" \M to t\ of the Polyakov line around a maximum 
of action density rises with increasing non-staticity St- 
In the region where one can separate the constituents ac- 
cording to the action density, the dipole moment is small, 
emphasizing again the homogeneity of the Polyakov line 
around the central point. In the region beyond S^ the 
dipole moment gradually stabilizes around a value of 1.5. 

In the same way as described so far, we have ana- 
lyzed configurations obtained by cooling at higher ac- 
tion plateaux. As an example we show in Fig. the 
histogram of non-staticity St for the same three temper- 
atures represented by N t = 4, N t — 5 and N t = 6. Al- 
though the precise border between static and non-static 
has not the clear meaning as for the one-caloron case, 
the trend is the same: at lower temperature the lumps of 
action tend to be more localized also in Euclidean time 



(" instanton-like" ) . Compared with the one-caloron case, 
the histograms are more smeared. 



So far we have investigated the outcome from cooling 
for varying temporal lattice extent and fixed spatial vol- 
ume. What happens when changing /3 at fixed lattice 
size and fixed lattice asymmetry ? In Figs. IHlwe again 
show non-staticity 5 t histograms without applying cuts 
for |L as |. We compare (a) for Nt = 4 the cases (3 = , 
f3 = 2.04 and (3 = 2.20 with each other, and (b) for 
N t = 6 the case (3 = 2.20 with f3 = 2.36 . [23 The 
statistics for the new /J-values (0.0, 2.04 and 2.36) is of 
order 0(200) field configurations each. Contrary to what 
one might have expected, we see that among the cases (a) 
and the cases (b), respectively, there are no qualitative 
differences. This means that the ratio of probabilities for 
the occurence of dissociated calorons (or separate static 
dyon pairs) vs. non-dissociated (single lump) calorons is 
not a function of the temperature describing the original 
equilibrium ensemble. Instead, it is mainly determined 
by the geometry (the aspect ratio) of the lattice. The 
larger the temporal lattice extent is in comparison with 
the spatial extent the lower is the probability to find sep- 
arated lumps of more or less static objects after cooling. 
In other words, once there is only one topological Q ^ 
object left, everything is classical. As a rule, the object 
is smooth, such that only the size of the box matters or 
other "infrared forces" (e.g. those forbidding single in- 
stantons) play a role. This observation clearly shows the 
limitation of the cooling method applied to the lowest 
action plateaux. 
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FIG. 6: The Polyakov line in the neighborhood of maxima of the action density : (a) modulus of the average summed-up 
Polyakov line \Ltot\ and (b) modulus of the corresponding "dipole moment" \M to t\ vs. non-staticity St, for the subsample with 
asymptotic holonomy near zero. For details of the definition see the text. The thick vertical line marks the non-staticity 8* 
where the caloron recombines. 
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FIG. 7: The same as in Fig. 3 (a) for S = 3Si ns t plateaux. 
No cut in the asymptotic holonomy has been applied. The 
thick vertical line marks the non-staticity 8* where the caloron 
recombines. 



IV. INSTANTONS OR CALORONS ON A 
SYMMETRIC 4-TORUS 



With low statistics (100 configurations ) we have also 
cooled equilibrium configurations generated with (3 = 2.2 
on symmetric lattices (16 4 representing "zero" tempera- 
ture). In this case we have found for the classical con- 
figurations at the plateau Sw ~ Sinst a broad distri- 
bution of non-staticity with a maximum around St w 2 



and with a tail extending beyond 3 (see Fig. |5J. These 
are obviously configurations with an action (topological 
charge) density well-localized in all four Euclidean direc- 
tions. There is a non-trivial behaviour of the Polyakov 
line inside these non-dissociated, instanton-like objects 
resembling the Polyakov line associated with the time di- 
rection in the finite temperature case for non-dissociated 
calorons. We have mapped the cooled lattice configu- 
ration with the help of all 4 possible definitions of the 
Polyakov line, which are now, on a symmetric lattice, 
physically equivalent to each other. Fig. EI] shows the 
profiles of action density, topological charge, and of the 
Polyakov lines (for four possible definitions) as they are 
seen in appropriate planes intersecting the lump through 
the centrum. The latter is defined as the maximum of 
the AD action density. For all types of Polyakov lines the 
characteristic double structure is seen exactly when the 
" asymptotic" value of the respective Polyakov line is not 
close to ±1. 

It is interesting to compare the pattern of the Polyakov 
line of these "caloron candidates" with the analytic KvB 
caloron formally constructed on a 16 4 lattice correspond- 
ing to maximally non-trivial holonomy with respect to 
what has been chosen as "time" direction. For this con- 
struction, the two constituents have been placed along 
the z-direction, separated by 8 lattice spacings. The 
lattice caloron is obtained calculating link by link from 
the continuum gauge field A^. Of course, such a con- 
structed lattice caloron has irregularities at the boundary 
if the lattice action is evaluated under the assumption of 
periodicity. After some cooling the configuration turns 
into an (approximate) solution of the lattice equations of 
motion on the torus. By that time the boundary arte- 



8 



N_t=4 
beta = 2.2 
beta = 2.04 
beta = 




0.2 



0.4 
deltaj 

(a) 



0.6 




FIG. 8: The distribution of non-staticity S t after cooling for fixed Nt but varying /3-values. No cut in the asymptotic holonomy 
is made. The spatial lattice size is always 16 3 . (a) N t = 4 : (3 = 0.0, 2.04, 2.20, (b) N t = 6 : (3 = 2.20, 2.36. 




FIG. 9: The distribution of non-staticity St after cooling for 
(3 = 2.2 on the symmetric lattice 16 4 , obtained from 42 con- 
figurations with unit topological charge. 



facts have disappeared (see Fig. Illjl . and the Wilson 
action has become Si ns t = 27r 2 /3. The KvB caloron is 
now adjusted to periodic boundary condition also in the 
x, y, z - directions. Despite the large separation the con- 
stituents formally have, judging according to the action 
density it is a non-dissociated caloron. It is seen that 
the cooling does not influence significantly the structure 
of the Polyakov lines of the KvB caloron witnessed by 
the time-directed Polyakov line (pit) showing the dou- 
ble peak structure mentioned above. The space-directed 
Polyakov lines (plx, ply, plz) have a simple structure char- 
acteristic for trivial holonomy. It is analytically clear that 



for caloron solutions asymptotically plx, ply, plz — » ±1. 
This is in contrast to the would-be "caloron candidate" 
obtained by cooling from confining equilibrium lattice 
configurations at zero temperature (Fig. I1UII where the 
double peak structure is present for all t, x,y,z - direc- 
tions and, generically, the asymptotic holonomy in all 
directions is non-trivial. 



V. CONCLUSIONS AND PERSPECTIVES 

In the present investigation we have subjected equilib- 
rium lattice gauge fields corresponding to various temper- 
atures to ordinary relaxation (usually called " cooling" ) in 
order to obtain an ensemble of classical solutions for fur- 
ther study. In this way we have extracted lowest-action 
classical solutions of unit topological charge typical for 
the given equilibrium ensemble. Notice that, at least as 
long this technique is used, this possibility is restricted 
to the confinement phase. We do not claim to find the 
real and complete topological structure hidden under the 
quantum fluctuations and consisting of topological lumps 
of compensating sign. What we wanted to see were the 
simplest solutions suitable as building blocks for a semi- 
classical modelling of the Yang-Mills path integral. In 
conclusion we can say, that cooling of equilibrium lattice 
fields in SU (2) lattice gauge theory shows that there ex- 
ist topological objects with a dyonic substructure, that 
we are able to resolve only in the confined phase. With 
increasing temporal lattice extent it becomes more proba- 
ble that the observed dyons recombine into calorons such 
that it becomes impossible to perceive the substructure 
looking exclusively at the distribution of action and topo- 
logical charge. Strictly speaking, we have found that the 
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FIG. 10: Profiles of the action density (s), the topological charge density (q), and of the Polyakov lines (pit, plx, ply , pi z) 
calculated along all straight line pathes parallel to the four axes for a 16 4 lattice caloron found by cooling a Monte Carlo 
generated equilibrium gauge field down to the one-instanton action plateau. The center of the caloron (at the maximum of its 
action density) was found at the site (a;, y, z, t) = (7, 8, 8, 14). The planes shown in the figures cross just this point. 
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FIG. 11: Profiles of the action density (s), the topological charge density (q), and of the Polyakov lines (plt,plx,ply,plz) as in 
Fig. 1101 but for a 16 4 lattice caloron obtained from a discretized KvB solution corrected by some cooling to the one-instanton 
action plateau (309 cooling steps). The center of the caloron (at the maximum of its action density) was placed to the site 
(x, y, z, t) = (8, 8, 8, 1). The planes shown in the figures cross this point. 
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fraction of dissociated calorons among all single-calorons 
events docs not depend on the physical temperature of 
the equilibrium fields but on the geometric aspect ratio of 
the box, i.e. the asymmetry of the lattice. However, all 
calorons have a nontrivial holonomy which is mapped out 
by the behavior of the Polyakov line inside and outside 
these configurations. 

In the limiting case of zero temperature {i.e. on a sym- 
metric lattice) topological lumps with Q = ±1 obtained 
by cooling look instanton-like and, at the same time, have 
the characteristic double peak structure of the Polyakov 
line in all t, x, y, z directions. This distinguishes them 
from the real KvB solutions which possess this structure 
only with respect to a distinguished "time" direction. 

Strictly speaking, an analytic solution of the Euclidean 
equations of motion with Q = ±1 is impossible on the 
AD torus 0. Nevertheless, on the lattice quasistable so- 
lutions of this kind exist. If considered only as lumps 
of action and topological charge, there is no contra- 
diction with the previous observations from "instanton 
searches" . As our analysis shows, for temperatures much 
lower than the deconfinement temperature rotationally 
symmetric (in AD) and (anti)selfdual lumps seem to be 
preferred under cooling. 

At finite temperature, they can be subsumed under the 
general class of KvB solutions. For zero temperature, 
however, a so far unknown parametrization with non- 
trivial asymptotic holonomy has yet to be found. 

After having completed the present investigation we 
were informed by Ch. Gattringer and R. Pullirsch [l6| 
about their paper "Topological lumps and Dirac zero 
modes in SU(3) lattice gauge theory on the torus" prior 
to publication. In this paper the authors concentrate on 
and more systematically continue the inspection of the 
low-lying modes of the chirally improved Dirac operator 
on the 4-torus, in the background of equilibrium lattice 
gauge fields at T = 0. It is extremely interesting that 
they find, for a certain fraction of Monte Carlo configu- 
rations in the \Q\ = 1 sector, a similar pattern of hopping 
zero-modes as a function of varying fermionic boundary 
conditions as for T ^ 0. Moreover, the change of lo- 
calization happens independently of which of the four 
directions is chosen as "imaginary time" direction along 
which periodicity can be purposedly changed. The con- 
current zero mode positions (two or three as in the finite 
temperature case) are consistent with being randomly 
distributed in the AD periodic box. The authors argue 
that this observation hints at the existence of a semiclas- 
sical background consisting of localized (in AD) instanton 
constituents. Whereas other (monopole ?) properties of 
the constituents are less obvious, a non-integer topolog- 
ical charge of the constituents has been hypothetically 
assumed in analogy with the dissociated KvB solutions. 
In fact, this should not be too difficult to be established. 



It is in particular this last interpretation that has to pass 
further tests. In case it becomes confirmed then it would 
be difficult to reconcile this with our observations based 
on cooling. Our findings are consistent with a picture in 
which (at finite temperature) the temporal size of the 
box determines (inversely) the size of the background 
solutions. In contrast, the scenario of Ref. |l(| seems 
to imply a complete dissolution (within the available 4- 
volume) of some constituents which still span a coherent 
semiclassical background. This would mean that there is 
no scale of coherence which is dynamically generated and 
decoupled from the overall size of the box. If the latter 
picture can be confirmed we would have to blame cooling 
(or the cooling action we have used) for artificially driv- 
ing all semiclassical configurations into integer-charged 
topological lumps at lower temperatures whereas sepa- 
rated monopole constituents are correctly reproduced by 
cooling only at higher temperatures. This would explain 
why previous studies of the topological properties using 
various cooling (smoothing) techniques (which actually 
missed the relevant temperature range) were not suit- 
able to discover the non-trivial holonomy accompanying 
all topological charges. 

There is still an ongoing debate on the question in as 
far a semi-classical interpretation of the QCD vacuum is 
valid at all 0}. For the time being we have nothing to 
add to this discussion. But our results point to the fact 
that the instanton gas or liquid model (as well as other 
models) have to be reconsidered by taking into account 
the non-trivial holonomy structure. We believe that this 
will improve the bad performance of (trivial holonomy) 
instantons in comparison with lattice smoothing results 
as reported in [l||[n|. 
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